Vortex Dynamics, Resistivity Formula, and 
Fluctuation-dissipation Theorems * 



X.-M. Zhu" and P. Ao^ 
Department of Experimental" /Theoretical'' Physics 
Umea University, S-901 87, Umea, SWEDEN 

Abstract 

We investigate the problem of forces on moving vortex in a superfluid or 
superconductor. The main purpose is to locate the source which leads to the 
contradictory results in the literature. We establish the connection between 
this problem to the difficult but well studied subject of resistivity formula 
in transport theory. The relaxation time approximation used in the force 
calculation via force-force correlation function is shown to be invalid. The 
roles of Berry's phase and fluctuation-dissipation theorems are discussed. 
For a vortex moving in a superfluid or superconductor, the only agreement we have 
reached so far is at zero temperature, in absence of any impurity potentials. Question arises 
in more complicated cases. Using the Berry's phase calculation [|l|], as well as the exact 
total force- force correlation calculation , It has been shown that the transverse force is the 
Magnus force, coming from extended states with no dependence on the core level spacing 
and other extrinsic details. Other microscopic derivations using core state transitions 
have shown, however, that the transverse force on a moving vortex is greatly reduced in 
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magnitude when the core level spacing is less than the energy scale associated with random 
scattering. Apparently, it seems that they are not two equivalent ways to obtain the same 
quantity but rather at least one of them must be incorrect or incomplete. Thus if we could 
settle the difference on the transverse force on moving vortex between core or extended states 
originated, we should be able to reach an agreement. However, these two ways of calculating 
the transverse force on moving vortex are shown to be equivalent to each other p]. It is 
impossible for them to give different results unless there is a hidden mistake elsewhere. 

The algebra in the derivation of vortex dynamics is either involved or abstract. The 
repeating of those calculations is not the most efficient route to detect and understand the 
source of disagreements. Instead, we note that the vortices satisfy a classical Langevin 
equation with parameters determined microscopically. This equation has the same form as 
a classical electron moving in a magnetic field. Therefore general properties of Langevin 
dynamics, for instance the fluctuation-dissipation theorems, should be obeyed by vortices. 
In addition, our knowledge on electronic transport can be readily used to help understanding 
vortex dynamics. In particular, we will examine the relaxation time approximation and show 
it to be invalid in resistivity calculation via force-force correlations. 

Let us consider a classical charged particle in a magnetic field obeying a generalized 
Langevin equation: 



Here the index i = x or y, the velocity of the particle u(t) = {ux(t) , Uy(t)) , m the mass, 
K(t) = {Kx{t), Ky{t)) external force, f(t) = {fx{t), fy{t)) a random force which simulates 
the effect of the thermal reservoir. The Einstein convention of the repeated indices as 
summation has been used. Beik Uk{t) represents the transverse force, the Lorentz force 
u(t) X B in the Langevin equation with the magnetic field taken along the 2;-direction. The 
matrix r]{t — t') = {f]ij} represents frictions in both longitudinal and transverse directions 



mui{t) 



- Jldt' 7]ik{t ~ t')Uk{t') + Ki{t) 

+Be,k Uk{t) + fi{t) . 
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of the particle motion. Its possible finite off-diagonal elements will change the effect of the 
original Lorentz force on the particle. In addition, we have 

< Mt) > = , 

< u^ito)fjito + t) > = 0, t > , (2) 



The second one is due to the causality and the last one is the equipartition theorem. 

The problem of particle response to a perturbation can be formulated into two different 
ways. We can calculate the velocity of the particle while the applied force is given. Here 
the Hamiltonian of the particle is known. In such a case, it is equivalent to obtaining a 
conductivity formula. Otherwise we can consider a given trajectory of particle motion and 
calculate what is the applied force needed to maintain it. It is equivalent to obtaining a 
resistivity formula, i.e. calculating electric field needed to maintain the given current. The 
derivation of vortex dynamics belongs to the second kind, where we assume a steady motion 
of vortex and calculate what is the external force acted on the vortex. In electron transport, 
we have a choice to formulate the problem in either way. In vortex dynamics, we can only 
formulate in terms of resistivity because the effective vortex Hamiltonian is unknown and is 
exactly what we want to obtain. 

The conductivity may be obtained by the Nakano-Kubo's formula, a calculation of 
velocity-velocity correlation function. It may also be obtained by solving Boltzmann equa- 



tion in presence of an electric field ||TT|. Both of them are standard methods. The relaxation 
time approximation is a valid one in these cases. 

The resistivity formula is much more confusing. Since 1960s, various resistivity formulae 
have been proposed and examined repeatedly |]9|-[12|. In the derivation of vortex dynamics, 
total force-force correlation formalism has been explicitly shown by Simanek to be the one 
used in various Green's function calculations of forces on moving vortex Q]. In addition, 
the force-balance theory has been used to calculate forces acting on a vortex p|. Recently, 



Boltzmann equation has also been used In the hterature of electron transport, all of the 
above methods have been used in obtaining a resistivity formula, although normally discussed 
without magnetic field. It has been demonstrated that the longitudinal resistivity calculated 
directly using total force-force correlation function is always zero in the zero frequency limit 
p!0| . Thus any formula that allows to obtain finite friction, or longitudinal resistivity, directly 



from total force- force correlation function must be incorrect in this limit. The longitudinal 
resistivity formula using Boltzmann equation gives either 0/0 results or a formula trivially 
equivalent to reciprocal conductivity formula ||lT . 



We will first briefly review conductivity formula. As expected, we will show that in 
such a case < u(t) > is determined by velocity-velocity correlation functions in the absence 
of K. We will verify that indeed the relaxation time approximation is valid in such a 
calculation. Parallelly, we will derive a resistivity formula by assuming a given velocity u{t) 
and calculating the average force < K > which should be applied to sustain such a motion. 
We find that the problems involved in vortex dynamics become clear after careful study of 
this simple model. 

Now let us put K = in Eq.(l) and calculate the velocity- velocity correlation function. 
Multiply Eq.(l) by Ujito) and take the ensemble average: 

m < Uiitf) + t)Uj{tQ) >= 

- /o dt' riikit - t') < Uk{h + t')Mi(to) > 

+Beik < Uk{to + t)uj{to) > + < fi{tQ + t)uj{to) > . 

The term < Ui(tQ)fj(to + t)> vanishes according to Eq.(2). 
Introducing a Laplace transform 

POO 

r^[co] = / dt e-'-'vit) , 
Jo 

defining the velocity-velocity correlation function matrix 

Uij{t) =< Uiito + t)uj{tQ) > , (3) 



and integrating by part, we have 

U[u] = {imuj + r][u\ - iayBy^ksT . (4) 

We have used the identity mU{0) — ksT. ay — —i{€ij} is the Pauh matrix. 
Next we calculate the total force-force correlation function matrix 

J^ij{t) < Ui{to + t)uj{to) > . 

Taking the Laplace transform, using the translational invariant in time 

< Ui{to + t)Uj{to) >= - < Ui{to + t)Uj{to) > , 

and the total force-velocity correlation function 

m < Ui{to + t)uj{to) > [u] — —mlAijiQ) + imujUij[u}\ , 

we have 

J-[uj\ — {iB(Ty -\- imuj -\- {mwYx 

{imw + r][uj] — iayB) j ksT . 

In the limit a; ^ 0, we have 

= iBuy UbT , (6) 

which is independent of r]. 

We then calculate the random force-force correlation matrix 

'R'ij{t)-<Mto + t)fj{to)> . (7) 

Prom Eq.(l) we can express TZ{t) in terms of total force-force, total force- velocity, and 
velocity-velocity correlation functions. Taking the Laplace transform and integrating by 
part, we obtain 
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n[uj] = r][Lu] mU{0) = r][uj] keT , (8) 

or ri(t) = IZif) I {ksT). This is the 'second' fluctuation- dissipation theorem described by 
Kubo 0. The generahzed friction ri{t) is given by the random force-force correlation. Be- 
cause the random force is determined by the thermal bath, 77 (t) has no off-diagonal part if 
the random force- force correlation matrix has not. An important conclusion we can draw 
from here is that the even though there is no time-reversal symmetry in the particle motion, 
the frictional force is always longitudinal as long as the heat bath does not generate an off- 
diagonal element in the random force-force correlation function matrix. For example, this is 
the case for a charged particle dynamics described by a single relaxation time in Boltzmann 
equation. 

Now we look for the connection between the correlation functions and the transport 
coefficients. First we consider the mobility. With an applied external force K(t) = K[co']e*'^* 
the mobility ii[uj\ is defined by < uju;] >= /ijj[ti;]Kj[u;]. From Eq.(l) we immediately obtain 
the mobility 



uluj] = {imuj + ri[uj\ — iayB) 
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in the limit to ~^ —00. Using Eq.(4), the mobility is related to the velocity- velocity correla- 
tion function 

= (9) 

This is the 'first' fiuctuation-dissipation theorem described by Kubo p, equivalent to the 
Nakano-Kubo's formula for the electrical conductivity. 

Next we consider that the particle is moving at a given velocity u(t) and find out what 
is the external force needed to sustain such a motion. It is equivalent to the calculation of 
resistivity if the particle is charged. From Eq.(l), we have the average force 

< >= {imuj + ri[uj] — iayB)ij Uj[uj] (10) 



which is trivially equivalent to the reciprocal of conductivity formula. Obviously this process 
does not provide us an independent way of calculating resistivity. 

However, if we are only interested in the average force < K > in a steady state motion, 
we do have an alternative resistivity formula. After taking cu ^ and using Eqs.(6) and (8), 
Eq.(lO) gives 



Taking to be a scalar (proportional to a unit matrix), the external force can take a more 
suggestive form. 



where the longitudinal component depends on TZ[0], the random force-force correlation func- 
tion, and the transverse component only on J-[0], the total force- force correlation function. 
Eq.(ll) is the steady state resistivity formula. It provides a direct way to obtain DC resis- 
tivity from force correlation functions. The straight forward interpretation of Eq.(12) is the 
force-balance: The externally applied force to keep the steady velocity is equal in magnitude 
but opposite in sign to the sum of the frictional and the Lorentz forces. Above analysis shows 
that the transverse force is not affected by the thermal reservoir under the assumption that 
77 [a;] is a scalar. 

So far, all of our calculations are exact. Now we will discuss how the results may change 
when employing the relaxation time approximation to account for the existence of thermal 
reservoir in comparison with the exact calculations. Without the thermal reservoir, the 
velocity-velocity correlation is given by 



Then we switch on the thermal reservoir to allow the relaxation process to happen. We 
use a relaxation time approximation by the standard rule, icu ^ icu + ri[Lj]/'m. Substituting 




(11) 




(12) 



L{[uj] — {imuj — iayB) U{0) . 



(13) 
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it into Eq.(13), we have found the velocity- velocity correlation under the relaxation time 
approximation is given by 

U[uj] — {imuj + r][uj\ — iayB)~^kBT , 

which is exactly the same as the one obtained by above rigorous calculation. We conclude 
that the relaxation time approximation can be a valid one for velocity-velocity correlations 
when used in a conductivity formula. 

Next we evaluate the force-force correlation by the relaxation time approximation. With- 
out thermal reservoir, the random force correlation is zero, that is, 7^(t) = 0. If we switch 
on the thermal reservoir by using a relaxation time approximation iuj —>■ iu) + r][Lu]/m, the 
random force correlation is still incorrectly set to zero, and cannot be made finite. The total 
force correlation without thermal reservoir is 

•^iiM — {iB(Ty -\- imuj -\- imujY 

j\. i \ y ^^^^ 

{imuj - iayB)-^)^f^mL{kjiO) ■ 
When we switch on the thermal reservoir by using a relaxation time approximation iuj — > 
iuj -\- r][uj\/m in Eq.(14), we have 

J-[uj] = {iBay + imu + ri[uj] — {imuj + r][u}]Y 
{imuj -\- r][uj\ — iayB)~^^ ksT . 

This is a rather complicated expression. We can simplify it in the limit uj « r/[a;], or 
UJT « 1, T = m/rjlO] is a relaxation time: 

•^t^] = m ^2 i^oT + iay{u;orf) ksT , (15) 

with iOf) = B/m. Let us use the resistivity formula Eq.(ll) to calculate the external force 
needed to keep the particle moving with a given velocity. With TZ[uj\ = Q and .^-'[0] given by 
Eq.(15), the external force is 

< K[0] >= ,2 ^0] + u[0] X B) . (16) 
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These results have no connection at all to the rigorous results shown in Eq.(12). Evidently 
the relaxation time approximation cannot be valid in such a calculation. 

However, Eq.(16) is familiar to us. With merely a re-definition of constants B = up with 
K the circulation and p the superfluid density, and c^o as the the core level spacing, this force 
becomes exactly the same as the one appeared in the derivation of vortex dynamics using 
the relaxation time approximation 0-0], where the average is made for all other degrees of 
freedoms except thOse of the vortex. The literal interpretation of such results is that the 
Lorentz force on a moving particle is canceled by the effect of thermal reservoir. It also shows 



a different friction |jT3[. The similarity of structure between Eq.(16) and those obtained in 
vortex dynamics with relaxation time approximation suggests that they may have the same 
source of error. The calculation by Simanek has explicitly used the relaxation time ap- 
proximation in the force-force correlation function. Although several other publications 
have not explicitly specified their methods as force-force correlation function calculations, 
their final formula are the same as Ref. [Q. 

There are also exceptions. In Ref. 0] , Stone made an attempt to use Boltzmann equation 
to solve the forces on a moving vortex. When we check Eq. (5.8) in Ref. 0, we find that if 
we substitute < k >=< k >o into this equation, where < k >o is the equilibrium value of 
< k > in the frame where vortex is stationary, the equation is not satisfied. However, < k > 
should relax to < k >o in such a case because there is no extra applied pinning force and the 
vortex is stationary in this frame. There are also additional problems. The equation itself is 
not Galilean invariant. Nonetheless, its solution has been transformed into lab frame using 
Galilean invariance in order to obtain force on a moving vortex . 

Before we conclude, we return to fiuctuation-dissipation theorems and Berry's phase. 
The type of fiuctuating forces we use in stochastic process will not generate any addition 
transverse force in low frequency limit because of their vanishing off-diagonal correlation. 
However, even though they are widely used, we may still question the validation of these 



9 



fluctuating forces. Our understanding is that wlien we separate a fluctuating force from what 
we leave to be systemic, we generally choose to assign the average effect to the later. Thus 
the fluctuating force no longer has a zero frequency component of off-diagonal correlation. 
The Berry's phase is more powerful in this respect. If indeed we have left an off-diagonal 
correlation with non-zero mean to 'random force', the zero frequency component will be in- 
cluded if we evaluate Berry's phase. Finally, we point out that when correctly evaluated, the 
core state transitions reproduces the results of Berry's phase calculation at zero temperature 



with impurities included, and the friction arises naturally fgJTj]. 
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